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Strong subadditivity inequality for a three-particle composite system is an important inequality 
in quantum information theory which can be studied via a four-particle entangled state. We use two 
three-level atoms in A configuration interacting with a two-mode cavity and the Raman adiabatic 
passage technique for the production of the four-particle entangled state. Using this four-particle 
entanglement, we study for the first time various aspects of the strong subadditivity inequality. 
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I. INTRODUCTION 

Entanglement 0] in a composite system refers to cer- 
tain implicit correlation between the subsystems arising 
from their interaction. It is the key resource of quan- 
tum computation and quantum information processing 
0- Due to recent advances in this field, entanglement 
has generated renewed interest. There have been dif- 
ferent approaches to understand and to quantify entan- 
glement Q- But so far the entanglement, only in a bi- 
partite pure state has been investigated very extensively. 
The von Neumann entropy 3 of either of the subsystems 
provides a good measure of entanglement in this case 
This is the quantum partner of the Shannon's entropy 
in classical information theory and is defined as 



Sii) = -Trj(pjlog2Pj) 



(1) 



where j = A, B. Here, pj is the reduced density operator 
of the subsystem j and is given by 



Pj = TripAB 



(2) 



where pAB is the density operator of the composite sys- 
tem under consideration and j, I ~ A, B, j ^ I. In gen- 
eral, the quantities S{j) satisfy the following inequality 
(due to Araki and Lieb) 0: 

\SiA) - SiB)\ < SiA, B) < S{A) + S{B) , (3) 

where S{A, B) is the joint entropy of the composite sys- 
tem comprising A and B. The second part of the above 
inequality is known as subadditivity inequality 9j . For a 
pure state, S{A,B) = and thus S{A) = S(b). The 
equality sign in the above relation holds good if and 
only if the composite density matrix pab can be written 
as a tensor product of its two reduced density matrices 
PA and pb, i-C, for a disentangled state. One can de- 
fine the index of correlation Ic given by the expression 
S{A) + S{B) - S{A,B) [13, which can also be inter- 
preted as information entropy in quantum information 
point of view. We note that Kim et al. have calculated 
the entropies of different kinds of pure states including 
two-mode Fock states and squeezed states JTj. Further, 
the above relation for entropy has been studied in the 
context of entangled Gaussian states [T^ . 



So far we have discussed about the measurement of 
entanglement in a bipartite pure state. If the compos- 
ite system is in a mixed state (defined by the density 
operator p), the entanglement of formation Ep can be 
defined in terms of the average von Neumann entropies 
of the pure states of the decompositions Q- Wootters 
has shown the quantity Ep to be an explicit function of 
p '^l- He has introduced the notion of concurrence in 
this context. 

We further notice that from the Schmidt decomposi- 
tion of a pure bipartite state, one can properly identify 
the entanglement present in the state This is also 

very useful to study bipartite continuous systems [T5l |. 
On the other hand, for a mixed state p, the separability 
criterion has been proposed ^1^] to study entanglement. 
This is based on positive partial transpose mapping of 
p. Thus the negativity (entanglement monotone) of the 
eigenvalues of the partial transpose of p could be a mea- 
sure of entanglement in a mixed bipartite system |17| . 
The concept of negativity as an entanglement measure 
has been used in context of interaction of atoms with 
thermal field The separability criterion has been 

extended to continuous systems ,1^ also. 

Despite many approaches to define entanglement for a 
bipartite system, there have been only a few approaches 
to quantify entanglement in the composite systems of 
three or more particles jTHMIP- 

We note that a gen- 
eralization of Schmidt decomposition in multipartite sys- 
tems in pure states has been introduced 22] . Coffman 
et al. 23| proposed a measurement of entanglement in 
a tripartite system in terms of concurrences of the pairs 
of subsystems. This measure is invariant under permu- 
tations of the subsystems. An average entanglement in a 
four-partite entangled state has been defined in terms of 
von Neumann entropies of the pairs of subsystems (2^ . 
Very recently, Yukalov has addressed the question more 
generally and quantified multipartite entanglement psl ] 
in terms of the ratio of norms of an entangling operator 
and of a disentangling operator in the relevant disentan- 
gled Hilbert space. 

In this paper, we put forward a possible measurement 
of entanglement of a four-particle system by studying 
the entropy of the reduced three-particle system. As 
mentioned above, the von Neumann entropy is a good 
measure for entanglement in a bipartite system. For a 



2 



tripartite composite state, this entropy satisfies a strong 
subadditivity inequality (SSI) j^j, which has many im- 
portant imphcations in the subject of quantum informa- 
tion theory. In this paper, we study the properties of a 
four-particle entangled state through the three-particle 
entropy and the SSI. 

The structure of the paper is as follows. In Sec. II, 
we provide a brief discussion on strong subadditivity in- 
equality from the quantum information point of view. 
In Sec. Ill, we describe a physical model and show the 
preparation of a four-particle entangled state. In Sec. IV, 
we study the validity of the SSI in the present context 
and provide a physical explanation of the results. We 
conclude this paper by proposing a measurement of the 
corresponding four-particle entanglement. 

II. STRONG SUBADDITIVITY INEQUALITY 

We have already mentioned that for a bipartite com- 
posite system of two particles A and B, the joint entropy 
S{A,B) satisfies the subadditivity inequality (PJ. For a 
composite system of three particles A, B, and C, this 
inequality can be extended to the following form (26j|: 

SiA,B,C) + S{B)<SiA,B) + SiB,C) . (4) 

This inequality is known as strong subadditivity inequal- 
ity. The most obvious situation that the equality sign 
holds in Q is when the composite density matrix pabc 
can be written as the tensor product of its three reduced 
density matrices as pA® PB ® pd i-e., when the system 
is in a disentangled state. However, the more stringent 
condition for this reads as 

l0g2 (pabc) - log2(pAs) = l0g2(psc) - log2(ps) ■ (5) 

There have been numerous implications of the above 
inequality I0J in quantum information theory Q . Firstly, 



it refers to the fact that the conditioning on the sub- 
system always reduces the entropy, i.e., S{A\B,C) < 
S{A\B), where, SiA\B) = S{A, B) - S{B) is the entropy 
of A conditional on knowing the state of B. Secondly, the 
above inequality implies that discarding a quantum sys- 
tem never increases mutual information, i.e., S{A : B) < 
S{A : B, C), where, SiA : B) = S{A) + S{B) - 5(A, B) 
is the mutual information of the subsystems A and B. 
Thirdly, quantum operations never increase mutual in- 
formation of two subsystems. This means that if the 
mutual information of the two subsystems A and B be- 
comes S'{A : B) after trace-preserving operation on B, 
then S'{A : B) < S{A : B). Further, this inequality 
Q implies that the conditional entropy of the subsys- 
tems A, B, and C is also subadditive, i.e., S{A,B\C) < 
S{A\C) + S{B\C). 

To verify SSI, one needs to calculate the entropies like 
S{A, B, C) which clearly requires a three-particle mixed 
state which we can produce using a pure four-particle 
entangled state Hj]. In the next section, we discuss how 
one can prepare a pure four-particle entangled state so 
that we can study SSI for the first time for a system 
realizable using cavity QED methods. 



III. PREPARATION OF FOUR-PARTICLE 
ENTANGLED STATE 

We consider two three-level atoms (A and B) with rel- 
evant energy levels in A-configuration (see Fig. ^ inter- 
acting with a two-mode high quality optical cavity. The 
specified annihilation operators for the cavity modes are 
a and b. The atoms are interacting with the cavity mode 
a in |e) <-*■ \g) transition and with the mode h in |e) <-> |/) 
transition. 

The Hamiltonian for the system under rotating wave 
approximation can be written as 



J 



H ^h{ujia)a + uj2b'^b) + h ^ [^eugAek){ek\+ ujfkgMk){fk\ + {gik\ek){gk\a + g2k\ek){fk\b + h.c.}] , (6) 

k=A,B 

where, toi^^mu is the atomic transition frequency between the levels \lk) and |mfc), ujj {j = 1,2) is the respective 
frequency of the cavity modes a and b, gjk (j = 1,2) provides the atom-cavity coupling. We assume gjk^s to be real 
and function of time. 

We start with the initial state \ipi) = \gA, gBT^, p), where n and p are the initial numbers of photons in the cavity 
modes a and 6, respectively and the two atoms are in \g) state. The state of the system can be expanded in terms of 
the relevant basis states in the following way: 

l'0(i)) = ci\gA,gB,n,p)+C2\gA,eB,n~l,p)+C3\gA,fB,n-l,p + l) 

+C4\eA, fB,n-2,p+l) + C5\eA,eB,n- 2,p) + C6\fA,fB,n- 2,p + 2) 

+C7\fA,eB,n- 2,p+l) +cs\fA,gB,n~ l,p+ 1) +Cg\eA,gB,n- l,p) . (7) 
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FIG. 1: Level diagram of two three-level atoms in A-configuration, interacting with two cavity modes defined by annihilation 
operators a and b. g^k and g2k (k — A,B) are the atom-cavity coupling terms for the k-th atom. A is the common one-photon 
detuning of the fields. 



From the Schrodinger equation we find the following equations of the corresponding probability amplitudes: 

di = -i{\fngiBd2 + V^giAdg) , 

d2 = -iiVngisdi + Ad2 + VM+l52BC?3 + Vn- IgiAd^) , 
ds = -f lg2Bd2 + Vn~ IgiAd^) , 

d4 = -i{Vn- IgiAda + Adi + ^/ + 1525^5 + ^/i + 252^1^6) , 

^5 = -i[Vn- l{giAd2 + gisda) + + '^{92Bdi + g2AC?7) + 2Ad5] , (8) 
= + 2g2Adi + ^/i + 2(/2bc?7) , 

= -f 132^^5 + + 2.925 + Ady + \/n - Igisds) , 

ds = - IgiB^T + + 1.92^1^9) , 

dg = -i{y/ngiAdi + Vn- Igisd^ + VM+l52yid8 + Adg) , 

where, we have used the following transformations: 

ci = di , 026"*-^* = d2 , C3 = da , €4^^* = d^ , cge^^'^* = ds , 

cg = de , cye"*"^* = d7 , cs = ds , cge"'^* = dg, cje^/fc - t^i,2 , (9) 

where, Zfe — gk, fk, is the one-photon detuning of the cavity modes for the fc-th atom. Here we have assumed that 
the cavity modes are in two-photon resonance and A^ = As = A. 

Writing these equations ^ in the matrix form [di] — — i[i\f][di], we find that one of the eigenvalues of the matrix 
[M] is zero. The corresponding eigenstate is 

iV'o) = [a\gA,9B,n,fJ,) + (3\fAjB,n-2,fi + 2) ~j\gA,fB,n- 1) - S\fA,gB,n~ 1,^+1))] , (10) 

I 



where 

a = (/2A.92S + + 2) , P = giAgiB\/n{n - 1) , 
7 = 9iB92A\/n{fi + 2) , S = giAg2B\/n{^j. + 2) , (11) 
P = ^Jo? -I- /32 -)- 72 ^ 52 ^ 

Clearly, this state is an entangled state of four parti- 
cles, namely, the atoms A and B, and the two modes a 
and 6. Using appropriate time-dependence of the pulses, 
the four-particle system can be prepared in this state, 
as discussed in the next section. Recently, there have a 
few experimental demonstrations of preparation of four- 
particle entangled state 28] and performance of a C-NOT 
gate I23. Interestingly, the state |^o) is a two-atom two- 



mode multipartite coherent population trapping (CPT) 
state which is a counterpart of the well-known CPT state 
for a single atom interacting with two coherent fields [s^l . 



IV. STUDY OF STRONG SUBADDITIVITY 
INEQUALITY 

We next discuss how the state |i/'o) can be prepared by 
using Raman adiabatic passage technique. We assume 
that both the atoms are initially in \g) state. We further 
assume the time-dependence of the Rabi frequencies gjk 
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of the two modes as 

9iA^9iB = gioexp(-(i-r)2/r^) ; 
52A=ff2B = g2oexp(-tVT^) • (12) 

Here, gjo [j — 1, 2) is the amphtude of the respective 
pulse, T and T are the width and time-separation respec- 
tively, of the two pulses. Note that the pulses are ap- 
plied in counterintuitive sequence. Under this condition, 
the atom-cavity system follows the evolution of the state 
l^o) adiabatically. This state is a zero eigenvalue eigen- 
state (adiabatic state) of the Hamiltonian 0. During 
this process, known as stimulated Raman adiabatic tech- 
nique (STIRAP) [sij . the atom-cavity system remains in 
this state for all times. In the present case, at the end of 
the evolution, the population of both the atoms are si- 
multaneously transferred to the state |/). However, if the 
atoms are not in one-photon resonance, i.e., if Ar ^ 0, 
then this transfer process is not complete. This happens 
because the system does not remain confined in the null 
adiabatic state |i/;o) for At ^ [sij . 

We now investigate the validity of SSI for any trio of 
quantum systems in the present process. We can express 
this inequality for any three particles, namely, atom A, 
atom B, and cavity mode a with number n of photons 
out of the four-particle system under consideration as 

E = S{A, B) + S{A, n) - S{A, B, n) - S{A) > . (13) 

Here, S defines the joint von Neumann entropy of the 
relevant subsystems [see Eq. This can be calculated 
from the state ^ by tracing over the other subsystems, 
e.g., 

S{A, B) = -TvAsiPAB log2 pab) , (14) 

where, pab is the reduced density matrix of the atoms A 
and B and is given by 

PAB^^ynA\^{t)){^{i)\) ■ (15) 

We show the time variation of E in Fig. |3 Clearly, E{t) 
never becomes negative during the evolution and thus 
the SSI lUni) holds for aU times. 

From Fig. 121 one clearly sees that for At = 0, in long 
time limit, E becomes zero. This means that the sub- 
systems (A, B, and the mode a with photon number n) 
become disentangled. This happens because of complete 
adiabatic transfer of population to the level |/) of both 
the atoms at long time limit. The entire process can be 
written as 

\gA,9B,n,p) — > \fA,fB,n-2,fi + 2) . (16) 

We have shown the time- variation of the coefficients a/P, 
13/ P, -i/P, and 5/P [see Eq. (dJ] in Fig. El This figure 
reveals the above evolution according to the state |V'o) 
under the action of the pulses H12|l . But for At = 60, 
since complete population transfer does not occur, the 
system remains entangled in the state at long time 
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FIG. 2: Time evolution of the parameter E for Ar = 
(solid curve) and Ar = 60 (dashed curve). This clearly shows 
that the strong subadditivity inequality remains valid in the 
present physical situation. The parameters chosen here are 
n = 2, M = 0, gjoT = 15 (j = 1, 2), T = 4r/3. 




t/x 

FIG. 3: Variation of the coefficients a/P, 13/ P, 7/P, and 5/P 
with time. The parameters chosen here are n = 2, ^ = 0, 
Qjor = 15, T = 4r/3, and Ar = 0. 

limit. This is clear from the dashed curve of Fig. 13 as 
the equality E' = no longer holds at this time limit. 

Thus we can recognize the expression E [see Eq. (|13ll ] 
as a measure of four-particle entanglement in the present 
process. Precisely, > 0, where the equality sign holds 
good for the disentangled states. An increase in value of 
E refers to increase in entanglement. Thus, during the 
evolution, the system gets more entangled for At = 
than for At — 60. However, at the end of the evolution, 
the entanglement persists for nonzero At. We must em- 
phasize here that, the present definition of entanglement 
measurement satisfies all the relevant criteria, namely, 
(a) it is semipositive, i.e., > 0, (b) E = for an disen- 
tangled state, and (c) the function E{t) is continuous in 
time domain. 

Using our four-particle entanglement, we can also 
study the inequality involving the entropies of two 
particles, say, atoms A and B. They remain strongly cor- 
related during the evolution, as Ic remains much larger 
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than zero. At the end of the evolution, Ic becomes zero 
for At — which means that the subsystems become 
uncorrelated. In other words, the entanglement between 
them vanishes. 

We note in passing that for a four-particle GHZ state 
defined by 

IV) = ^(10000) + 11111)) (17) 

of four qubits A, B, C, and D, one is led to a three-particle 
mixed state pabc defined by 

PABC = i(|000)(000| + |111)(111|) . (18) 

In this case, S{A, B, C) = S{A) = S{A, B) = 5(5, C) = 
log2 2 — 1. Therefore, the parameter E in this case be- 
comes zero, as from Eq. (|13|l . So we have a counter- 
example, in which the equality sign in |0J) holds for an 



entangled state, too. However, we note that the above 
state (|18f) satisfies the condition (|3J) and thus the said 
equality. 



V. CONCLUSIONS 

In conclusion, we have shown for the first time the 
role of strong subadditivity inequality for entropies in a 
four-particle composite system. The stimulated Raman 
adiabatic technique has been used to prepare the four- 
particle entangled state using two three-level atoms ini- 
tially in their ground states in a two-mode cavity. We 
further show that the parameter E could serve as a pos- 
sible measurement of entanglement in the four-particle 
entangled state under consideration. 
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